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THE PRINCIPAL AXES OF STELLAR MOTION 
By H. Raymond 

DUDLEY OBSERVATORY. ALBANY. N.Y. 
Communicated by E. B. Frost. May 28. 1917 

Suppose the stars of known motion transferred to a fixed point O, 
and allowed to move with their present velocities for a definite length of 
time. They will have expanded into a cluster, the 'velocity-figure.' 
According to the older idea that the stars were moving 'at random,' 
this cluster should be spherical. Since J. C. Kapteyn showed"^ the 
figure to be elongated, several methods of finding the amount and 
direction of elongation have been proposed and used. Some of these 
make use only of the directions of observed proper-motions, leaving 
the amounts of those motions wholly out of account, thus avoiding the 
difficulties caused by our ignorance of the distances of most stars; and 
are based on some hypothesis as to the actual shape of the velocity- 
figure. The most notable of these are A. S. Eddington's^ and K. 
Schwarzschild's.* The method here used makes no assumption as 
to the form of the velocity-figure, and is essentially a process of finding 
its principal momental axes. It thus belongs to the same genus as 
those used by C. V. L. CharUer,* K. W. GyUenberg,^ A. S. Eddington 
and W. E. Hartley,* and H. C. Plummer.' 

Let the coordinates of a star in the velocity-figure, referred to suitable 
axes through 0, — say to the First Point of Aries, to (6h, 0°), and to the 
North Pole, respectively, — be x' , y' , z'; the direction cosines of an 
arbitrarily chosen line be I, m, n; and the projection of the star S upon 
this line be P. Then the component of motion in the direction I, m, n, 
is p = OP = Ix' + my' + nz'. Representing the mean of a quantity 
by enclosing it in square brackets, 

[pp] = P [x'xl+m^ [y'y']+n' [z'z'] +2lm [x'y'] +2 In [x'z'j +2mn [y'z'] (A) 

The axis of preference is defined as the line for which [pp] is a maximum; 
for the axis of avoidance [pp] is a minimum. 

Newcomb, in his paper 'on the Principal Planes toward which the 
Stars Tend to Crowd,' applies to another problem an essentially similar 
process, and gives a detailed solution by Lagrange's method for a 
maximum or minimum of (.4), from which the procedure of this investi- 
gation is adapted. Putting A = [x'x'], B = [y'y'], C = [z'z'], D = 
[x'y'], E = [x'z'], F - [y'z'], the maximum, minimax, and minimum 
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values of {pp\ are the three roots of the cubic 
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D {B-\) F 

E F (C - X) 
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Wx']^[xx]-[x]-\x] 


[x'y'] = [xy] - 


-W-M 


[y'y>] = [yy] - \y] . [y] 


[x'z'] = [xz] - 


-M-W 


[z'z'] = [zz] - [z] . [z] 


[y'A = bz] - 


-M-W 



and the corresponding values of /, m, n, are derived by inserting these 
roots successively in the set of equations 

(.4-X)/+ Dm+ En = 

Dl + {B-\)m + i?w = (C) 

EI+ Fm+{C-\)n==0 

The coordinates x', y', z' are referred to the unknown center of gravity 
of the velocity-figure. We can observe only the velocities x, y, z, 
relative to the moving sun. But we readily find that 



(D) 



The quantities [x], [y], [z], are the negatives of the coordinates of solar 
motion, according to the method of Bravais. 

We are limited to the use of the proper-motions alone, or the radial 
velocities alone, since to combine them requires a knowledge of indi- 
vidual parallaxes which we do not possess. The rectangular components 
of the proper-motion part (tangential to the celestial sphere) of the 
stellar motions, resolved as before, may be expressed in terms of the 
proper-motions m, h', and the spherical coordinates of the stars, while 
IX, IX may be expressed in terms of x, y, z. Designating the tangential 
part by subscript I, we readily obtain 

Xi = x{l —cos^a cos^8)—y sin a cos a cos^6 — z cos a sin a cos 5, 

and similarly for yi and Zi. From these x, Xi, etc., are found by multi- 
pUcation. The means of these quantities are found by integration over 
the sphere. They are, cleared of fractions, 

3M = 2[x] 15[xiyi] = 7[xy] 

3W = 2b] 15[a;iZi] =7[xz] 

3W=2[z] lS[yiZx] = 7bz] 

15[a;ia;i] = 8M-f- [yy] + [zz] 

15[yiyi]= [xx] + 8[yy]+ [zz] 

15 [zizi] = [xx] + [yy] + 8[zz] 



(E) 



The admissibility of these integrations depends upon certain con- 
ditions, which may be stated in the form of three h)T)otheses: 
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I. The stars are evenly distributed over the sky. 

II. The distribution in distance is the same in different parts of the 
sky. 

III. The motions of the stars are distributed in the same way in all 
accessible parts of space. 

These hypotheses are of course not rigorously true. Ill is open to 
objection, but if it is not at least approximately true we have no general 
problem, merely a number of local problems. After the mean square 
velocities in the three principal directions are known, failure of II, so 
far as effect of different mean square parallax in different parts of the 
sky is concerned, may be taken care of and a new approximation made. 
The theoretical value of [xiXi] + [yiyi] + [ziZi] can be computed for each 
region of the sky and compared with the observed values. A difference 
may be considered as due to the stars being nearer or farther than 
the average and factors applied accordingly. These approximations 
converge very rapidly; in fact, the error from the source considered is 
small, even in the first solution. Errors due to the spread or range 
of the stars in distance are not corrected by this process, but must 
usually be negligible. Hypothesis I can be made valid by simply 
weighting the material proportionally to the area of the sky covered, 
instead of proportionally to the number of stars. 

The material used was the 5943 proper-motions less than 80" per 
century from Boss's Preliminary General Catalogue. Of these, 5384 
were less than 20" per century (allowing for solar motion*) and were 
called "S"; the remaining 559 stars of larger proper-motion were called 
"L." Group S was subdivided according to type of spectrum as fol- 
lows: B, types Oe5 to B 5, 492 stairs; A, B 8 to A 4, 1647; F, A 5 through 
F, 656; G, 446 stars; K, 1227; M, 223; X, including stars of unknown 
type — mostly faint — and some of types N, 0, and P, 693 stars. 

The numerical results appear in the accompanying tables. Table 1 
gives, in column 1, the designation of the group; in 2 and 3 the Right 
Ascension and Declination respectively of the direction in which [pp] 
is a maximum — the vertex of preferential motion; in 4 and 5, the same 
for direction of least motion. The third direction, in which [pp] is a 
minimax, is necessarily at right angles to each of these, hence is not given. 
The last three columns give the value of [pp] for these three directions. 
Unit of motion, 1" per century. 

Table 2 gives the solar motion derived at the same time. A and D 
are the Right Ascension and Declination of the Apex, M the amount 
of solar motion. The last three columns are given as a means of com- 
paring the various groups. If the means without regard to sign of the 
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components of motion in the three principal directions be [u], [v], [w], 
we have 

[u] = V (2Xi/t), [v] = V (2X2/7r), M = V (2X3/X) 

provided the components are distributed according to the error law. 
They have merely the broad characters of such distribution, — sufficient 
to justify the process in the same sense that the customary method of 
finding the 'probable error,' without inquiring too dihgently as to 
whether the errors of observation are actually distributed in the manner 
assumed in theory, is justified. 

TABLE 1 

Systematic Motion 





POLE OF PKEFERENCE 


POLE OF AVOD>ANCE 


MEAN SQUARE MOTIONS 




Al 


Di 


A3 


Dj 


Al 


X2 


X, 


B 


100?6 

92.4 
93.1 
93.8 

93.1 
94.0 
83.9 

91.9 
92.5 
92.3 


-38?4 

+ 1.8 
+2.1 
+3.4 

+ 11.3 
+25.3 
+30.2 

+ 4.2 
+35.3 
+ 16.9 


180?4 

182.9 
185.2 
185.4 

191.2 
205.1 
192.6 

183.8 
211.7 
184.7 


+ 12?6 

+ 15.2 
+45.9 
+24.4 

+35.3 
+37.2 
+29.0 

+24.9 
+34.5 
+ 8.1 


4.58 

25.65 
,58.05 
32.26 

29.30 
28.14 
21.08 

28.91 
745.6 
107.3 


1.94 

9.15 
28.78 
15.81 

18.71 
12.58 
12.72 

14.99 
261.4 

58.7 


1.33 


A 


2.57 


F 


22,66 


G 


9.96 


K 

M 


11.46 
3.48 


X 


5.33 


S 


9.89 


L 


194.7 


All 


31.8 







TABLE 2 

Solar Motion 



GROUP 


NUMBER 


SOLAR APEX 


MEAN MOTIONS 




A 


D 


M 


(ul/M 


(vI/M 


(wl/M 


B 


492 

1647 
656 
446 

1227 
223 
693 

5384 

559 

5943 


274?2 

266.7 
262.9 
256.6 

272.5 
272.2 
271.4 

269.0 
269.3 
269.1 


+39?3 
+27.9 
+27.0 
+45.8 
+38.4 
+45.7 
+42.1 

+32.1 
+31.5 
+32.3 


2 "33 
4.03 
6.36 
2.60 
4.02 
3.45 
3.06 

3.83 

22.17 

5.42 


0.73 
1.00 
0.96 
1.74 
1.08 
1.23 
1.20 

1.12 
1.01 
1.53 


0.47 
0.60 
0.67 
1.22 
0.86 
0.82 
0.93 

0.81 
0.61 
1.13 


0.40 


A 

F 


0.32 
0.60 


G 


0.97 


K 


0.67 


k 


0.43 


X 


0.60 


s 


0.66 


L 


0.52 


All 


0.83 
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In regard to these results the following points seem especially worth 
noting: 

1. All groups show markedly unequal motions in three principal 
directions. In all groups the direction of greatest motion, or axis of 
preference, lies near the plane of the galaxy not far from its intersection 
with the equator, the axis of avoidance is nearly perpendicular to the 
galaxy. 

Using radial velocities Gyllenberg* finds much the same axis of pref- 
erence, but the other two axes practically interchanged. So do Edding- 
ton and Hartley* for later tj^es, but their results for early t3^es agree 
better with table 1. Are these differences real, in the sense of being 
tj^ical of what may be expected whenever radial velocities and proper- 
motions are compared? If so, what peculiarity of distribution or of 
motions in space, out of the many that might easily be suggested, lies 
back of it? Here is a subject well worth investigation. 

2. The separation of the poles of preference into two groups corre- 
sponding respectively to 'early' and 'late' types,' is confirmed. There 
seems to be some evidence of progression of the same character within 
each of the two groups of types, as if with advance of tj^e the 6-hour 
vertex moved northward. With advance of type there is an increase 
of stellar as compared with solar motion; that in the preferred direction, 
however, increasing less than the other two components. 

3. The apices of solar motion show the division into two groups al- 
ready found by several investigators; A and F having apices in smaller 
R. A. and Decl. than K, M, and X. B seems to resemble the latter 
group. G is anomalous. 

These results agree rather closely with those found by L. Boss by 
Airy's method,** and by the writer by Schwarzschild's method,' all 
three being based upon practically the same material. This agreement 
will be found to extend even to details. 

4. The distribution of stars of type B is so peculiar that to weight 
the material according to area covered would have been unfair. The 
integrations which give equations (E) can be made over zones parallel 
to the galaxy, and the resulting relations used to correct each zone for 
perspective separately, the material then being combined with weights 
proportional to the number of stars. The results in the tables were 
thus found; a check solution by the ordinary method gave for the vertex 
(122°,-36°). 

The vertex of preferential motion is near the antapex of solar motion. 
This may in part be a manifestation of the tendency, noted above, for 
the earlier types to have more southerly vertices. In part it is due to 
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the fact that the stars' 'spread' in distance tends to produce a false 
elongation of the velocity-figure in the direction apex-antapex. This 
tendency is stronger the greater the spread, the smaller the velocity- 
figure relative to solar motion, and the more nearly spherical the figure. 
The last two causes, and possibly the first, operate strongly in the B 
stars. There is no appearance of such an effect in the other groups. 

5. The separation of groups "S" and "L" was made some years ago, 
for a different purpose, and according to a method which made allow- 
ance for solar motion relative to group S, but not enough for L, and 
made none for preferential motion. Consequently, while S no doubt 
represents fairly well the smaller motions, L does not truly represent the 
large ones. In group "All," while the large motion characters are 
diluted by large numbers of small motion stars, the former have in- 
dividually much more weight than the latter; this therefore is probably 
more representative that L of the large proper-motion stars. A re- 
division into small and large proper-motions, on principles better adapted 
to the purpose of this investigation, is in progress. 

• London, Rep. Brit Ass., 1905, (257). 

2 London Mon. Not. R. Astr. Soc, 67, 1906, (34). 
^Gottingen, Nachr., Math.-Phys. KL, 1907, (614); 1908 (191). 

* Medd. Lunds Astr. Obs., (Set. 2), No. 9, 1913. 
^Ibid., No. 59, and (Ser. 2), No. 13, 1914-15. 

<> London Mon. Not. R. Astr. Soc, 7S, 1915, {521). 
Ubid., 76, 1915, (121). 
Ms/r. /. Boston, 612, 1910, (98). 
UUd.,616, 1915. 
^"Ibid., 614, 623-4, 1910-11. 



